arXiv:1509.01439v2 [physics.optics] 25 Sep 2015 


Controlling circular polarization of light emitted by quantum dots using chiral 

photonic crystal slab 

S. V. Lobanov,i’2’3 s. G. Tikhodeev,!’^’^ N. A. Gippms,^’! A. A. Maksimov,^ E. V. Filatov,^ 

1. 1. Tartakovskii,^ V. D. Kulakovskii,^ T. Weiss,^ C. Schneider,^ J. GeBler,^ M. Kamp,^ and S. Hofling^’^ 

^A.M. Prokhorov General Physies Institute, Russian Aeademy of Seienees, Vavilova Street 38, Moseow 119991, Russia 
^Sehool of Physies and Astronomy, Cardiff University, Cardiff CF24 3AA, United Kingdom 
^Skolkovo Institute of Seienee and Teehnology, Novaya Street 100, Skolkovo 143025, Russia 
^Institute of Solid State Physies, Russian Aeademy of Seienee, Chernogolovka 142432, Russia 
^M. V. Lomonosov Moseow State University, Leninskie Gory 1, Moseow 119991, Russia 
^4^^ Physies Institute and Researeh Center SCoPE, 

University of Stuttgart, Stuttgart D-70550, Germany 

^ Teehnisehe Physik, Physikalisehes Institut and Wilhelm Conrad Rontgen Researeh Center for Complex Material Systems, 

Universitdt Wurzburg, D-97074 Wurzburg, Germany 

^SUPA, Sehool of Physies and Astronomy, University of St Andrews, St Andrews, KY16 OSS, United Kingdom 

(Dated: 25 September, 2015) 

We study the polarization properties of light emitted by quantum dots that are embedded in 
chiral photonic crystal structures made of achiral planar GaAs waveguides. A modihcation of the 
electromagnetic mode structure due to the chiral grating fabricated by partial etching of the wave¬ 
guide layer has been shown to result in a high circular polarization degree pc of the quantum dot 
emission in the absence of external magnetic held. The physical nature of the phenomenon can be 
understood in terms of the reciprocity principle taking into account the structural symmetry. At the 
resonance wavelength, the magnitude of \pc \ is predicted to exceed 98%. The experimentally achieved 
value of \pc\ — 81% is smaller, which is due to the contribution of unpolarized light scattered by 
grating defects, thus breaking its periodicity. The achieved polarization degree estimated removing 
the unpolarized nonresonant background from the emission spectra can be estimated to be as high 
as 96%, close to the theoretical prediction. 

PACS numbers: 71.36.+C, 42.65.Pc, 42.55.Sa 


I. INTRODUCTION 

The possibility to control the polarization state of ra¬ 
diation from quantum emitters has drawn attention of 
researchers in recent years as it opens various important 
applications in spin-optoelectronics, quantum informa¬ 
tion technology, chiral synthesis and sensing. A common 
way for polarization conversion and rotation is the use 
of wave plates of birefringent materials and optical grat¬ 
ings. An imbalance between the left- and right-circularly 
polarized photons takes place in chiral materials with 
nonequivalent left- and right-circularly polarized electro¬ 
magnetic field modes .1^® The imbalance in semiconduc¬ 
tors is usually reached by applying a static magnetic field 
that disturbs the time-reversal symmetry and leads to a 
splitting of the left- and right-circularly polarized modes. 

Advances in nanoscale fabrication have made it pos¬ 
sible to realize artificial photonic structures with de¬ 
sired symmetry and density of the environmentally al¬ 
lowed electromagnetic modes, which allows to control the 
spontaneous emission rate of emitted light, its radiation 
pattern and direction. Moreover, the fabrication of chi¬ 
ral nanostructures from achiral semiconductor materials 
opens the possibility to control the polarization of emit¬ 
ted light ,1^ as the polarization effects are determined by 
the overall structure symmetry.!^ This method has consid¬ 
erable advantages in comparison with the others as it is 
compatible with semiconductor technology and allows to 


fabricate nanoscale devices (thickness of a conventional 
quarter-wave plate is much larger) with a simple opera¬ 
tion (external magnetic field application is not easy). 

Recently Konishi et al.^^ fabricated a simple chiral 
nanostructure that consisted of InAs quantum dots 
(QDs) embedded in the waveguide region of a planar 
GaAs based dielectric waveguide coupled to a dielectric 
chiral photonic crystal slab (GPCS), and obtained a de¬ 
gree of circular polarization of QD emission pc = 26%. 
A significantly higher value of pc of QD emission (up to 
70% in the direction normal to the structure plane and 
up to 81% at some angle) was experimentally achieved 
from a planar semiconductor microcavity with chiral par¬ 
tially etched top mirror,!^ but this structure is far more 
complex. The degree of circular polarization of emission 
from a chiral structure made of achiral semiconductors 
depends strongly on its geometry and the radiation fre¬ 
quency, which allows one to optimize the structures. Re¬ 
cent numerical simulations^ show that values of pc as 
high as 98 and 99% in waveguide and micro cavity struc¬ 
tures, respectively, can be reached by their optimization. 

In this paper, we have experimentally investigated the 
polarization of emission of InAs QDs embedded in a pla¬ 
nar GaAs based dielectric waveguide coupled to a di¬ 
electric chiral photonic crystal slab with sub-wavelength 
period in order to reach the predicted high values of pc- 
The structure geometry has been optimized on the base 
of the computation of the QDs radiation performed using 
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an oscillating point dipole modelP^^ A circular polariza¬ 
tion degree of In As QD emission in the finite (NA less 
than 0 . 02 ) aperture as high as 81% is experimentally re¬ 
alized in fabricated optimized structures. The achieved 
value is found to be limited by the unpolarized scattered 
light background caused by the fluctuations in the peri¬ 
odic potential due to the size and shape imperfections of 
etched rectangular pillars in the chiral layer. It is worth¬ 
while to note that the magnitude of pc obtained from the 
spectra with canceled background due to the lattice im¬ 
perfection reaches 96%. Thus, the considered very simple 
chiral structure made of conventional achiral semicon¬ 
ductor materials using the chiral morphology effect can 
enable the fabrication of compact solid state circularly 
polarized light-emitting devices with future progress in 
fabrication technology. 
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FIG. 1. (Color online) (a) Schematics of a structure with 
GaAs/air chiral photonic crystal slab, GaAs waveguide layer 
with an embedded layer of QDs, and a Alo. 7 Gao. 3 As cladding 
layer on GaAs substrate; (b) top view of one period of the 
right (or counterclockwise) twisted structure. 

In addition, we clarify the physical nature of the phe¬ 
nomenon. Based on the reciprocity principle and the 
symmetry of the structure, we show the possibility to en¬ 
hance and filter out one circularly polarized light compo¬ 
nent from the originally unpolarized QD emission. This 
is a consequence of the Fabry-Perot resonances between 
propagating photon modes in the chiral photonic crystal 
slab. We show that the circular polarization efficiency of 
such a chiral filter has to oscillate with its thickness, and 
the CPCS in the optimized structure is actually a half¬ 
wave plate. Its thickness can be made however much less 
than in the case of using natural optically active mate¬ 


rials with a weak chirality, of the order of the emitted 
light wavelength A only. It happens because in a dielec¬ 
tric CPCS with a high contrast of dielectric constants it 
becomes possible to achieve a large difference between 
wavenumbers of different propagating modes, of the or¬ 
der of = 27r/A. 

We would like to add that recently, a possibility to 
control deterministically the direction of photon emission 
(depending on the handedness of its circular polarization) 
by a QD inside a chiral photonic crystal waveguide has 
been demonstrated.^^ Our structure with randomly dis¬ 
tributed QDs may open an alternative way to this goal, 
being not sensitive to the exact positioning of the emit¬ 
ting QDs inside the structure. 

The paper is organized as follows: The experimental 
details and the sample are described in Sec. [H] . The ex¬ 
perimental results are described in Sec. |III| T he physical 
nature of the effect is discussed in Sec. IIVI Section El 
contains the concluding remarks. 


II. EXPERIMENT 

The epitaxial layer structure, which formed the basis of 
this study, consisted of a 691 nm thick GaAs membrane 
with a single layer of embedded randomly distributed 
In As quantum dots (Fig. The membrane is grown 
on top of a 1 pm thick AlGaAs buffer layer. We fabri¬ 
cated an array of nano-pillars into the upper nominally 
428 nm thick part of the membrane layer via electron 
beam lithography and reactive ion etching. The nano¬ 
pillars were arranged into a square lattice of unit cells 
consisting of four rectangular pillars, each rotated by 90° 
with respect to its nearest neighbors. The lattice period 
p was varied in the range of 0.7-1.1 pm. The unit cell 
is shown in Fig. Et- The distance between pillars was 
p/5, their sides were p/5 and 2p/5, respectively. Both 
right and left (mirrored) twisted pillars were manufac¬ 
tured. (For definite ness , we call hereinafter the struc¬ 
tures shown in Figs. | 1 | 2 ] as right-twisted.) The thickness 
of GaAs waveguide layer (ecaAs = 12.42), containing a 
plane with QDs was 263 nm. The plane with QDs was lo¬ 
cated 100 nm above its bottom. The lower Alo. 7 Gao. 3 As 
(^AlGaAs = 9.66) buffer layer was 1 pm thick. The whole 
system was grown by molecular beam epitaxy on a GaAs 
substrate. Representative SEM pictures of the right- 
twisted chiral photonic layer are shown in Fig. 

The sample was at a temperature of 5 K in the insert of 
an optical cryostat. The luminescence polarization was 
analyzed by a quarter wave retarder and linear polariz¬ 
ers. The angle resolution of 0.5° is ensured by a small 
aperture. The circular polarization degree of the lumi¬ 
nescence from the structure prior to etching was less than 
0 . 1 %. 
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FIG. 2. Representative scanning electron microscopy (SEM) 
pictures of the right-twisted chiral photonic crystal fabricated 
from the GaAs-based waveguide, top (a) and side (b) view. 


III. EXPERIMENTAL RESULTS 


Experimental circularly polarized spectra of QD emis¬ 
sion are shown in Fig. (left axis) for left and right- 
twisted rectangular pillars. The green lines in Fig. 
right axes, show the resulting spectra of the circular po¬ 
larization degree of emission 


I+ +1-' 


( 1 ) 


where is the emission intensity in right/left circular 
cr^-polarization.C^The spectra are recorded in the direc¬ 
tion 2 ; (perpendicular to the structure plane) with nu¬ 
merical aperture NA=0.005. For comparison, the cal¬ 
culated emission spectra in polarization and pc are 
shown in Fig. |^for the right-twisted structure with the 
same parameters. The spectra were calculated using the 
reciprocity principle and the optical scattering matrix 
method.l^ Compared to Ref. [H we changed the param¬ 
eters of the materials, because the current experiment is 
carried out at low rather than at room temperature. 

Figure shows that, in accordance with the calcula¬ 
tions, the polarization of the emission line at the reso¬ 
nance frequency is opposite in the structures with the 
left and right-twisted pillars. 

The experimental and calculated dependencies in 
Figs.andare in a good qualitative agreement. A small 
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FIG. 3. (Golor online) Experimental circularly polarized 
spectra of QD emission from the fabricated chiral waveguide 
structures with (a) the right- and (b) the left-twisted rect¬ 
angular pillars. Green dashed lines show the corresponding 
circular polarization degree pc. 


9 nm red shift of the calculated resonance wavelength is 
most probably due to a slightly incorrect dielectric sus¬ 
ceptibilities used in the calculations, the agreement can 
be easily improved by their fitting. 


Note, however, that (i) the experimental line FWHM is 
about 5 times greater than the calculated one (^ 4 and 
0.8 nm, respectively), and (ii) the line sits on a back¬ 
ground whose magnitude is about 9% from its peak in¬ 
tensity whereas the background in the calculated spec¬ 
tra does not exceed 1%. Both the line broadening and 
background magnitude are found to depend strongly on 
the pillar quality, which indicates small deviations from 
the ideal structure. Figure shows that magnitudes of 
pc = —0.81 and -1-0.67 in the waveguides with the right 
and left-twisted pillars are limited by the unpolarized 
background. It is worthwhile to note that the magnitude 
of I Pc I at the resonance reaches 0.96 when this unpolar¬ 
ized background is removed as shown in Fig. This 
value is very close to the calculated one, which indicates 
that the predicted values of |pc| ~ 0-98 could be reached 
due to further advances in fabrication technology. 
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IV. DISCUSSION 


A. Angular dependencies of emission intensities, 
Pc^ and the quasiguided modes 


As already discussed}^ the peaks in the QD emission 
spectrum are due to the excitation of lossy quasiguided 
modes in the GaAs waveguide inside a periodically mod¬ 
ulated structure. Their expected angular dispersion is 
shown in Fig. where false colors show the calculated 
intensity of the emission (in (t~ polarization) from the 
right-twisted structure with period p = 770 nm as a func¬ 
tion of energy (wavelength) and polar angle ^ in M — F 
(diagonal) and T — X directions. These Bragg resonances 
can be understoocP^ as corresponding approximately to 
the folding of the dispersion of the guided modes into the 
1st Brillouin zone (BZ) of the 2D grating (square lattice 
in our case), 


= —J{k, + IKbz? + {ky + mK^zY, ( 2 ) 


where k = {kx^ky) is the in-plane wavevector, K^z = 
27r/p is the Bragg wavenumber, p is the period, Ueff is 
the effective refractive index of the guided mode, I and 
m are integers. It appears that the resonance of interest 
at Ao ^ 901 nm in the current structure corresponds to a 
fourfold degenerate /, m = ±2 state in the F-point (in the 
scalar approximation and neglecting the energy splitting 
due to the point symmetry group), see Fig. It is also 
in agreement with the analysis of Ref. |9l see, e.g., the 
calculated in-plane distribution of the emission intensity 
in the dominant polarization in Fig. 6 b of Ref. |9l 

The degeneracy of the state ^± 2,±2 is lifted at /c 7 ^ 0. 
According to Eq. near the F-point 


EprnO^) — Eq -\- 


hc{lkx + mky) 




( 3 ) 



FIG. 4. (Color online) Calculated circularly polarized spec¬ 
tra and circular polarization degree of QD emission from the 
chiral waveguide structure with the right-twisted pillars. The 
sign of circular polarization is opposite for the left-twisted 
pillars. The emission intensity is normalized to the emission 
intensity of the same oscillating dipole in vacuum (without 
the structure). 



FIG. 5. (Color online) Experimental circularly polarized 
spectrum of QD emission from the fabricated chiral wave¬ 
guide structures with the right-twisted pillars. Fig. 1^, with 
subtracted unpolarized background. 



Polar angle 0 in M-F-X directions, deg 


FIG. 6. (Color online) Dispersion of quasiguided modes in 
spatially structured GaAs waveguide. False color scale shows 
the calculated QD emission intensities (in left polarization, 
logarithmic scale) at various k = kosinO, ko = Ct;/c, along 
rX and FM (diagonal) directions in the 1st BZ. Labels /,m 
(see in the text) mark the order of the corresponding Bragg 
resonance in the F-point. 


where Eq is the resonance energy in the F-point, Eq = 
E± 2 ,± 2 {k = 0). Thus, in the 0 — A, 0 — T directions, ac¬ 
cording to Eq.[^ the (± 2 , ± 2 ) modes split into two degen¬ 
erate doublets with nearly linear /c-dispersion E — Eq = 
Ehckx^y/{V^nQf[). In the diagonal directions, the modes 
with I = m = ±2 split into two singlets with nearly lin¬ 
ear dispersion E — Eq = EhckjnQ^^ whereas those with 
I = —m = ±2 is a degenerate parabolic doublet with 
an effective mass rueff ^ 10 “^me, where rrie is the free 
electron mass. Thus, the modes with I = —m = ±2 are 
nearly dispersionless at small k <0.2 As a result, 

the lines of equal energy at E = Eq are two mutually 
perpendicular lines at ky = Ekx in this range of k. That 
is well seen in Eigsj^,d representing the measured an- 
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FIG. 7. (Color online) Polar plots of the measured (left panels) and calculated (right panels) angle dependencies of the QD 
emission intensities in left and right circular polarizations. The measured spectra were recorded for the resonance wavelength 
Ao = 901.1 nm (panels c and d), corresponding to the (±2, ±2) Bragg resonance a,t E = Eq (see in the text). Panels (a,b) 
and (e,f) show the corresponding dependencies for Ao — 3.4 and +2.6 nm. The simulated spectra were calculated around the 
theoretical resonance wavelength Ao,th = 910.5 nm (panels i and j). Panels (g,h) and (k,l) show the corresponding dependencies 
for Ao,th + 3 nm. 


gular distribution of the QD emission near ^ = 0 in two 
circular polarizations at the resonance wavelength. 

Figures [^,b and e,f show that these lines at 
split into doublets ky = ^ Ak and ky = —kx^ Ak 

in accordance with the D 4 point group symmetry of the 
modes’ dispersion at the F point. 

The measured anglular distributions of the emission in¬ 
tensities in polarizations and of the circular polariza¬ 
tion degree (left panels of Figs.[^and [^respectively) are 
in a good qualitative agreement with the calculated ones 
displayed in the right panels. It is seen in Fig. that, 
in agreement with the calculation, the strongest emis¬ 
sion in the dominated (cr“) polarization is at the cross¬ 
ing of the isoenergetic lines and reaches its maximum at 
E = Eq when the state at /c = 0 becomes fourfold de¬ 
generate. According to calculations shown in Fig. the 
QD emission intensity in the dominating left polariza¬ 
tion at the resonance E = Eq exceeds that in vacuum 
more than 30 times, whereas that in the opposite polar¬ 
ization is strongly (about 3 times) suppressed. The effect 
is connected to a modification of the environmentally al¬ 
lowed electromagnetic modes by the CPCS in /c-space. 
As a consequence, it is very sensitive to the quality of 
the photonic crystal. Any disturbance of the potential 


periodicity results in the increased Rayleigh scattering 
contributing to /c = 0 emission in the suppressed polar¬ 
ization. The result is the decreased degree of circular 
polarization in the experiment. 


B. Symmetry analysis and the physical origin of 
circularly polarized emission 

Although the possibility to obtain nearly 100% cir¬ 
cularly polarized emission from unpolarized emission of 
randomly distributed QD was already demonstrated nu¬ 
merically,!^ the physical mechanism of this effect has not 
been yet clarified. The key property of the structure that 
helps to understand this mechanism is the expected the¬ 
oretically oscillatory dependence of the degree of circular 
polarization of emission of such structures on the thick¬ 
ness of the chiral photonic slab thickness, illustrated in 
Fig.H While the calculated dependencies of the emis¬ 
sion intensities in right and left polarization on L (see 
Fig. IT) consist of series of at first sight incoherent peaks, 
the dependence pc{L) (see green line in Fig. [^) looks as 
coherent Fabry-Perot oscillations. In what follows, based 
on the reciprocity and symmetry properties of our struc- 
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FIG. 8. (Color online) Polar plots of the measured (left 
panels) and calculated (right panels) angle dependencies of 
\pc\ of the QDs emission, corresponding to the data in Fig.[^ 


ture, we show that the chiral photonic crystal slab in the 
structure of interest plays the role of a wave plate. It 
introduces a phase shift, proportional to AKzL {L is the 
CPCS thickness), between the propagating modes with 
a large (of the order of ko = uj/c) difference of wavenum¬ 
bers AKz. We will show that these modes are basically 
different waveguide modes in the GaAs rectangular verti¬ 
cal waveguides, i.e., in the building blocks of the CPCS. 
There is also a large constant phase shift between left 
and right circular polarizations, due to the absence of 
horizontal mirror symmetry of the photonic crystal slab 
with different top and bottom boundaries (air/CPCS and 
CPCS/GaAs, respectively). As a result, the transmission 
of the CPCS in left and right polarizations oscillate as a 
function of L in counter-phase, which opens a way to 
reach nearly 100% polarization selectivity. 


1. Reciprocity principle for the structures with Ca and Da 
symmetry 

Th e em ission intensity of QDs can be calculated di- 
rectl^ffSES® or using the electromagnetic reciprocity 
principleIn the latter case, the emission problem of 
randomly distributed and randomly polarized oscillating 
dipoles is replaced by an illumination problem, in which 
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FIG. 9. (Golor online) The calculated dependencies of 
(a) the emission intensities and (b, green line) the cir¬ 
cular polarization degree of emission pc at the resonance 
wavelength Ao,th = 910.5 nm on the chiral photonic crys¬ 
tal slab thickness L. Blue line in panel (b) shows the cal¬ 
culated circular selectivity of the chiral slab transmission 
Pc,T = (T + - T-)/(T+ -h T ), see the discussion in Sub¬ 
sec. [I^ 


the structure is illuminated by a polarized plane wave 
Ein(^, y, z) = Ee^ ex.-p{ikoz — iut). (Hereinafter index j 
marks the polarization state, e.g., j = x^y in case of lin¬ 
ear polarization along x, respectively, and j = E for 
right and left polarization.) Then, the resulting electric 
field distribution in the layer with QD, 

P = lldxdy{\Eif + \Eif), (4) 

gives the corresponding polarized component of the QD 
radiation intensity.!^ 

The reciprocity principle helps to understand easily, 
why in the structure with a C 4 symmetry, the emission of 
QD in the vertical direction cannot be linearly polarized, 
and why for a structure with higher symmetry (e.g., D 4 ), 
it cannot be also circularly polarized. 

Assume E(x, y, z) and H(x, y, z) are the solutions of 
Maxwell’s equations for the system with C 4 symme¬ 
try. Then, the electromagnetic field rotated by 90° is 
also a solution. The 90° rotation switches the linear 
polarization of the illuminated plane wave j = x to 
the orthogonal polarization j = y, while the integral 
Q does not change, . Because of the Eq. 

Pxy = — I^)/EI^) = 0, hence the emission cannot 

be linearly polarized. 

Similarly, if a structure has a D 4 symmetry (equiva¬ 
lent to C 4 ^), Maxwell’s equations are invariant also under 
mirror reflection operation, e.g., in XOZ plane, which 
switches the sense of circular polarization of the illumi¬ 
nated plane wave f ± to the opposite one, j 
while the integral Q does not change, /+ = I~. Thus, 
again, the degree of circular polarization of emission 
pc = (/^ — /“)/(/+ = 0 , hence the emission cannot 

be circularly polarized, too. 
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2. S-matrix approach and Ca symmetry 

In the structure of interest with C 4 symmetry, a layer 
with randomly distributed and randomly polarized QD 
is placed in the planar GaAs waveguide. Two lower effec¬ 
tive dielectric permittivity layers operate as mirrors. The 
bottom mirror is the Alo. 7 Gao. 3 As layer. The top mir¬ 
ror is the 2D chiral photonic crystal. The propagation 
of light in this system can be described by the optical 
scattering matrix. The basic idea is (i) to split the sys¬ 
tem into several layers in which the dielectric permittivity 
does not depend on z, and (ii) to find the 2D electromag¬ 
netic modes of each slab, their electromagnetic field dis- 
tributions {Ex{x,y),Ey{x,y)}, {Hx{x,y),Hy{x,y)}, and 
wavenumbers Kz. The periodicity of the structure along 
X and y axes leads the fields to be represented as products 
of periodic functions and harmonics exp(i(/ca,x + kyp)), 
where kx^y are in the 1st BZ, with —i^Bz/2 < kx,y < 
KbzI‘^- 

The G 4 symmetry leads the in-plane components of the 
electric and magnetic fields to be represented after a 90° 
rotation in one of the four forms 

{E:,{-y,x),Fy{-y,x)} = {Fy{x,y),-E:,{x,y)} , (5) 

where n = —1,0,1,2, and F means the electric or 
magnetic fields. The advantages of the periodicity 
and C 4 symmetry can be simultaneously used only 
for kx = ky = 0. In this case, the left polarized light 
_ ^ex±^^ikoz-iujt ^ propagating in air normal to 

the structure plane, corresponds to n = — 1 , while the 
right polarized one corresponds to n = -hi. Since two 
other forms are optically inactive (see, e.g., in Ref. [T3|) . 
we will consider later only n = ± 1 . 


3. Time-reversal symmetry 

The time-reversal symmetry of Maxwell’s equations for 
systems consisted of non-absorbing and non-gyrotropic 
materials relates wavenumbers and electromagnetic fields 
of the propagating (Imi^;^ = 0 ) modes with oppo¬ 
site n. If z) and UsDix^y, z) are the solu¬ 

tions of Maxwell’s equations for infinitely thick layer, 
the time-reversal symmetry requires E 3 ^(x, z) and 
—H 3 ^(x,^, z) to be the solutions, too. 

Considering the propagation of some mode with fields 
{Ex{x,y),Ey{x,y)}e^^-^ and {Hx{x,y), Hv{x, y)} , 

one obtains that {E'* (x, ^), £’*(x, ^)} e and 

— {H*{x,y), Hy{x,y)} are also solutions. It 

corresponds to the counter propagating mode with op¬ 
posite polarization n, since complex conjugation of one 
polarization (n = -hi or — 1 ) switches the polarization 
to the opposite one (n = —1 or + 1 , respectively). As a 
result, 

Kt = = {H-)\ (6) 


where the upper indices ± show the values of n. Note 
that the eigenmodes are always defined up to an arbitrary 
factor C. If propagating eigenmodes are normalized to 
have fixed energy flow (we choose ^ per elementary unit 
cell), the eigenmodes are still defined up to an arbitrary 
factor with unitary absolute value e*^. Hereinafter, we 
choose the phase ^ for the opposite circular polarizations 
so that their electromagnetic fields are complex conjugate 
of each other. 


4- Unitary scattering matrix 


The optical scattering matrix can be approximately re¬ 
placed by the unitary scattering matrix S'^^.C^The thicker 
the chiral photonic crystal layer is, the better is this ap¬ 
proximation. The schematics of the QD emission near 
the Ao,th = 910.5 nm resonance with n = -|-1 or —1 (de¬ 
coupled in case of C 4 symmetry) is shown in Fig. 10 r. 
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FIG. 10. (a) Scheme of light scattering in the struc¬ 

ture. (b) Schematics of k = 0 light emission near the 
A = 910.5 nm resonance. Only the propagating modes 
are shown, with wavenumbers (in units of ko) Kzi" = 1; 
Kzi',2',3'a' = 2.61,1.85,1.32,1.22; K^i = 1.11. Actually, 
there are 12 more propagating modes in the waveguide layer 
(including the /c = 0, Kzjko = yT = 3.52 mode), but in com¬ 
parison with the quasiguided mode 1 they do not scatter back 
from the bottom interface of the waveguide even if they are 
excited on the top interface and thus can be neglected. 


Here, St and Sb are the unitary scattering matrixes 
of the air-GPGS and CPGS-waveguide interfaces, respec¬ 
tively. Hereinafter we omit the subscript u in the scat¬ 
tering matrix notation for simplicity. 

In general, the unitary scattering matrix expresses the 
amplitudes of outgoing modes in terms of the amplitudes 
of incoming modes: 


/ outi \ 


/ini\ 


/ ri ti2 ti3 . \ 

/ ini 

out 2 

— S' 

in2 


hi ^2 hs • I 

in2 

outs 

— 0 

ins 


hi h2 ^3 • I 

in3 

V ••• y 




V ../ 

V... 


( 7 ) 

The diagonal elements of the unitary scattering matrix 
are the reflection amplitudes of the channel, and 








































- Calculation with 401 harmonics 

Model with calculated papameters 
- Model with approximate parameters 



aaaaaa 

A 


0 2 4 6 8 10 



L, |am 


FIG. 11. (Color online) The calculated dependencies of the 
CPCS transmission in (a) right and (b) left polarizations 
on CPCS thickness L. Panel (c) shows the calculated circu¬ 
lar selectivity of transmission pc,T- Blue lines are the results 
of calculations with 401 spatial harmonics. Green and red 
lines correspond to the two-modes model Eq.[^with approx¬ 
imated and calculated (using 401 harmonics) parameters (see 
in Tab.|^, respectively. 


non-diagonal elements are the transition amplitudes tij 
from the channel to the channel. 

Similar to the relations in Eq. the time-reversal sym¬ 
metry of Maxwell’s equations for systems consisting from 
non-absorbing and non-gyrotropic materials relates also 
reflection and transition amplitudes in polarizations: 

4 = r-,t±=t-^. (8) 


5. Two-mode resonances 


In the structure of interest, at the resonance wave¬ 
length, the QDs excite mode # 1 , which leaks along 2 : 
axis primarily through the channels 7 ^ 1 ' and ^2' with 
transition amplitudes and respectively, see 


Fig. Wj. For both polarizations cr+ or cr“, these modes 
propagate through the CPCS with wavenumbers Kzv 
and Kz 2 '^ respectively. At the top boundary of the CPCS 
these modes transform into the mode 7 ^ 1 " with the tran¬ 
sition amplitudes and , respectively. The 

total transition probability, including the phase shift due 
to propagation over the CPCS, is 








1 ^ 2 ' 




(9) 

where L is the thickness of the CPCS. To simplify the 
formulas, let’s introduce the new variables 


= 4|t 


1 ^ 1 ' 


, /-T 


( 10 ) 


^ 1 ^ 2 ' 


^ 2 '^V‘ 






( 11 ) 


‘P = arg 


^1^2' ^2'^V‘ 




AK = Kz 2 ' -Kzv. Then 


0 

~ 4 


1 + T] e 


±iAKL+^^ 


If T]^ = 1, then 


± _ T’i „„„2 


cos 


( 12 ) 


(13) 


(14) 


is the oscillating function of the thickness L, which varies 
from 0 to T^. Moreover, if \cp~^ — (p~ = tt, the transmis¬ 
sion in one polarization reaches its maximum when the 
transmission in the opposite polarization reaches mini¬ 
mum and vice versa. The circular polarization selectivity 
of transmission pc,T = — T“)/(T+ + T~) in these 

thicknesses appears to be around 100 %. 

A simple model to estimate the introduced above pa¬ 
rameters for the resonance at Eq is given in the Ap¬ 
pendix. The results of the simplified model are compared 
with the optical scattering matrix calculation in Fig. pTj 
Although the parameters of the simplified model seem 
to differ substantially from the ones calculated within 
the scattering matrix method with 401 harmonics (see 
in Tab. [^, the results from this crude model for the cir¬ 
cular polarization selectivity of transmission (panel c in 
Fig. are in a good quantitative agreement with the 
calculations within the scattering matrix. It happens be¬ 
cause 77 ^, |(/ 9 + — ip~ I and especially the pc are sensitive to 
the ratios of the partial transmissions trough modes V 
and 2 ', which deviate less from the calculated values. 

The dependence Pc^t{L) is also shown as blue line in 
panel b of Fig.[^ It can be seen that it is very close to the 
dependence pc{L) (green line in the same panel). This 
fact supports our assumption that it is the selectivity of 
transmission through the chiral photonic slab that con¬ 
trols the degree of circular polarization of QDs emission 
in our structure. 


V. CONCLUSION 

In conclusion, we have fabricated and studied a struc¬ 
ture with light emitting InAs QDs inside a GaAs waveg¬ 
uide with a chiral photonic crystal layer on top that 
demonstrates a strongly circularly polarized photoemis¬ 
sion of QDs at zero magnetic field with \pc\ up to 81% and 
up to 96% if the unpolarized background due to CPCS 
imperfection is removed. It indicates that the predicted 
values of \pc\ close to 100% can be reached with further 
progress in fabrication technology. 

The measured angular dependence of the circularly po¬ 
larized emission intensities near the resonance is shown 
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to agree with the calculations using the optical scatter¬ 
ing matrix. This allows to attribute the circularly polar¬ 
ized resonant emission to the excitation by QDs of the 
(± 2 , ± 2 ) Bragg resonance of quasiguided modes in the 
GaAs waveguide. Based on the reciprocity and symme¬ 
try analysis of the structure, we show that the chiral pho¬ 
tonic crystal slab in the suggested waveguide structure 
works as a half-wave plate, exploring the Fabry-Perot in¬ 
terference between the propagating modes in CPCS, to 
reach nearly a 100 % circular polarization of the trans¬ 
mission. The advantage of using the photonic crystal 
slab with a large contrast of dielectric permittivities is 
its giant chirality, which allows to fabricate a very thin 
half-wave plate, with a thickness of the order of the emit¬ 
ted light wavelength. Additional advantage of the CPCS 
wave plate is its in-plane rotational isotropy, due to the 
C 4 symmetry. 
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Appendix A: Analytical estimation of transition 
amplitudes 


Here we investigate the resonance shown in Eig. 
in the optimized structure and assume the free space 
wavelength of light to be fixed and equal to 910.5 nm. 
The transition amplitudes ti^j in Eqs. ( 10 )-( 12 ) can_be 
roughly estimated via the normalized (see after Eq. 
in-plane components of the electric and magnetic fields 
of modes i^j as 




-IL 


'P/2 E* X Hj + Ei X H* 
da;d2/-^-^ 

-p/2 2 


(Al) 


Here, the symbol x means the pseudoscalar product 


Exli = E,Hy-EyH,. (A 2 ) 

The CPCS consists of dielectric rectangular pillars that 
are known to support the waveguide modes (in the verti¬ 
cal direction). The electric (magnetic) field distribution 
of the waveguide mode with the larger/smaller value of 
wavenumber can be roughly approximated by a ho¬ 
mogeneous field parallel to the longest/shortest (short¬ 
est/longest) side inside the pillar and zero field out¬ 
side it. Considering four rectangular pillars each ro¬ 
tated by 90° with respect to its nearest neighbors re¬ 
quires the electric and magnetic fields to satisfy Eq. (§. 


Note, due to the condition (§, the integral in Eq. |A 1 


can be taken only over the first quarter of the elemen¬ 
tary unit cell (it should be also multiplied by 4), since 
E*l-y^x) X = ^*{x,y) X B.{x,y). 

The electric and magnetic field distributions of modes 
# 1 ' and # 2 ' in the first quarter of the elementary unit 


cell are 



E^(a:,y) = j 

f TiCis 

if xe (|,f)&ye ( 0 ,f) 


otherwise. 



(A3) 

Upx,y) = 

lo> 

ifxe (|,f)&ye ( 0 ,f) 
otherwise. 


(A4) 


E2/(a:,2/) = 


^ r Ci£-V4e„ if x€ {I f)ky€ ( 0 , f) 


0, 


otherwise, 


(A5) 


h: 


= ( 2 ; y) ^ if a; e (|, ^) & y e (O, 


lo> 


otherwise. 


(A 6 ) 


where the constant Ci is introduced for brevity. 


2y/2p' 


The electric and magnetic field distributions of the 
modes 7^1 (which is the (± 2 , ± 2 ) Bragg resonance of the 
GaAs layer guided mode, in the horizontal direction) and 
# 1 " are 


T 

~7r~ 


Ef{x,y) = 


47r 


ko 


py K 


cos ( —{x ^ y) ) ± i 


. T 


^/2 




COS ( ^{x — y) 


P 


(A7) 


V 2 


Uf{x,y) = -J^ 
P\ ko 

47r A -e.^ + 

—{x + y) ±^ 


V2 


4:71 

-{x-y) 

p 

(A 8 ) 



(A9) 

(AlO) 

where the wavenumber is 


A'„ - pS ’p' . 

(All) 

5 = 12.42 is GaAs dielectric permittivity. 

Substituting (A3)-(A10) in (|A 1 |), one can 
sition amplitudes 

find the tran- 

iTi. =C2(-C3±iC4), 

(A 12 ) 

=<^2(-C4±iC3), 

(A13) 

y-± _ ^ 

(A14) 
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where the constants C 2 , C 3 , C 4 , and C 5 are 


^2 = 


IOtt^ 


£i/4j^+e-i/4 




zl 


K,i 

ko 


^ • 2 
C 3 = Sin 

5 

^ 27r TT 

64 = cos ^ + cos 
5 5 

C5 = i(ei/" + e-^/"). 


(A15) 

(A16) 

(A17) 

(A18) 


Substituting formulas ( |A12|)-([Al3 into the definitions 
of T(^, 77 ^, and i.e., Eqs?^ 0 ^ 11 ), and ( 12 ), respec¬ 
tively, one finds 


T^ = cici{ci + ci)^ 0.007, 


77^ = 1, 


C 3 


\ip~^ — if I = TT — 4atan— ^ 111°. 

64 


(A19) 

(A 20 ) 

(A 21 ) 


Thus, as discussed above in Eq. transmission in each 
polarization depends periodically on the CPCS thickness. 
As the function of thickness, it varies from 0 to 0.007 and 
is shifted by 111 ° with respect to the transmission in 
the opposite polarization. The comparison of this simple 
model with the numerical simulation with 401 harmonics 
is given in the Tab. |Tj The resulting dependencies of the 
CPCS transmission in cr+ and a~ polarizations and of the 


circular polarization degree of transmission in compari¬ 
son with the same formulas (but using exact parameters 
from Tab. and the results of simulation with 401 har¬ 
monics are given in Fig. It can be seen that this very 
approximate and crude model is in a qualitative agree¬ 
ment with the numerical simulations. The dependencies 
in Fig. calculated with 401 harmonics, demonstrate 
also sharp peaks, which are due to the influence of other 
electromagnetic modes in the CPCS, not taken into ac¬ 
count in this simple two-mode model. More detail on the 
comparison of the model with numerical simulations are 
given in the Supplementary material. 


TABLE I. Comparison of the simple model in Appendix with 
the numerical simulation with 401 harmonics 


Variable 

Approximate 

Simulation 

Relative difference 


-0.0256 + 0.0830i 0.0093 + 0.0673i 

25% 

IV' 

^ 1 ^ 2 ' 

-0.0256 - 0.0830i -0.0487 - 0.0564i 

22 % 

-0.0830 + 0.0256i -0.0452 + 0.0220i 

28% 

C^2' 

-0.0830 - 0.0256i -0.0417 + 0.0199i 

46% 

t + 

0.4820 

0.5774 - 0.0613i 

11 % 

2'^l" 

0.4820 

0.5734 - 0.0994i 

13% 

0.4820 

0.7290 - 0.0561i 

21 % 


0.4820 

0.7277 - 0.0122i 

20 % 


0.007 

0.0062 

6 % 

TV 

0.007 

0.0075 

4% 

77 + 

1 

0.932 

4% 

77 “ 

1 

0.775 

13% 

|(^+ - 

111 ° 

139° 

11 % 


^ Here, we calculate the relative difference S of two variables a 
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The spatial dependences of the scalar product in Eq. (Al) for the approximate 
, ^ 1 ^ 2 ' ’ ^ 1 ^ 2 '’ described in Eqs. (A3-A8), are shown in Figs. OIS, 038 
, 058 , 078, respectively. Figures 028, 048, 068, 088 show the corresponding 
dependencies, calculated within the optical scattering matrix with 401 harmonics. 
Left (right) panels show the real (imaginary) parts of the scalar product (Al). 

Figures 098, 108, 118, 128 show the in-plane spatial distributions of the electric 
and magnetic fields of the modes #1, #1', and #2'. 8ee the description of these 
modes in 8ubsec.IV.B.5 Two-mode resonances and Appendix A. 
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Figure OIS. Spatial dependence of the scalar product in Eq. (Al) for the 
approximate ■ Approximate analytic electric and magnetic fields from 
Eqs. (A3), (A4), {Al), and (A8). 



Figure 02S. Spatial dependence of the scalar product in Eq. (Al) for the 
approximate . Calculated electric and magnetic fields. 





















Figure 03S. Spatial dependence of the scalar product in Eq. (Al) for the 
approximate . Approximate analytic electric and magnetic fields from 
Eqs. (A3), (A4), (A7), and (A8). 
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Figure 04S. Spatial dependence of the scalar product in Eq. (Al) for the 
approximate . Calculated electric and magnetic fields. 
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Figure 058. Spatial dependence of the scalar product in Eq. (Al) for the 
approximate t7 . Approximate analytic electric and magnetic fields from 


Eqs. (A5)-(A8). 
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Figure 068. Spatial dependence of the scalar product in Eq. (Al) for the 
approximate • Calculated electric and magnetic fields. 




















Figure 07S. Spatial dependence of the scalar product in Eq. (Al) for the 
approximate ■ Approximate analytic electric and magnetic fields from 
Eqs. (A5)-(A8). 



Figure OSS. Spatial dependence of the scalar product in Eq. (Al) for the 
approximate ■ Calculated electric and magnetic fields. 
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Figure 09S. Electric (top panels) and magnetic (bottom panels) field distributions 
of the mode #1 in n = +1 (left panels, see also Fig. lOS) and n = — 1 (right 
panels) polarizations. See the explanation in Subsec.IV.B.5 Two-mode resonances 
and Appendix A. In all panels, the blue (red) arrows mark the real (imaginary) part. 
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Figure lOS. Real (left panels) and imaginary (right panels) parts of electric (top 
panels) and magnetic (bottom panels) field distributions of the mode #1 in n = +1 
polarization. See the explanation in Subsec.IV.B.5 Two-mode resonances and 
Appendix A. See also left panels in Fig. 09S. 
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Figure IIS. Electric (top panels) and magnetic (bottom panels) field distributions 
of the mode #1' in n = +1 (left panels) and n = —1 (right panels) polarizations. 
See the explanation in Subsec.IV.B.5 Two-mode resonances and Appendix A. On 
all panels, blue vectors mark real part and red vectors mark imaginary part. 
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Figure 12S. Electric (top panels) and magnetic (bottom panels) field distributions 
of the mode #2' in n = +1 (left panels) and n = — 1 (right panels) polarizations. 
See the explanation in Subsec.IV.B.5 Two-mode resonances and Appendix A. On 
all panels, blue vectors mark real part and red vectors mark imaginary part. 











































































































